A Causal Source which Mimics Inflation 
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How unique are the inflationary predictions for the cosmic microwave anisotropy pattern? In this 
paper, it is asked whether an arbitrary causal source for perturbations in the standard hot big bang 
could effectively mimic the predictions of the simplest inflationary models. A surprisingly simple 
example of a 'scaling' causal source is found to closely reproduce the inflationary predictions. This 
Letter extends the work of a previous paper (ref. 6) to a full computation of the anisotropy pattern, 
including the Sachs Wolfe integral. I speculate on the possible physics behind such a source. 

The prospect of mapping the cosmic microwave back- 
ground (CMB) anisotropies to high resolution has raised 
the exciting possibility of confirming fundamental theo- 
ries of the origin of structure in the universe. The in- 
flationary theory, being the simplest and most complete, 
is the present front-runner, and the latest CMB mea- 
surements do even seem to hint at the presence of first 
Doppler peak in the angular power spectrum, at I ~ 200 
as predicted by the simplest, spatially flat inflationary 
models. These spectra are distinct from those predicted 
by cosmic defect or baryon isocurvature models, and it 
is an important question whether spectra of this form 
are really a unique prediction of inflation. Or could a 
non-inflationary mechanism somehow replicate them? 

The fundamental difference between inflationary and 
non-inflationary mechanisms of structure formation is 
that inflation alters the causal structure of the early uni- 
verse, adding on a prior epoch during which correlations 
are established on scales much larger than the Hubble ra- 
dius. This is of course how the standard 'horizon puzzle' 
is solved. The perturbations produced during inflation 
are susceptible to very detailed test, and one could hope 
that the same feature of 'super-horizon' scale correlations 
(quotes indicate a standard big bang definition) could be 
used as a signature. If 'super-horizon' perturbations were 
shown to exist, this would strongly confirm inflationary 
structure formation, for no causal mechanism could have 
produced them within the standard big bang. 

Since COBE observed large scale perturbations on the 
CMB sky, one might think this already showed that 
'super-horizon' perturbations were present at last scat- 
tering. But these large angle anisotropies could have 
been produced causally within the standard big bang, 
through a time dependent gravitational potential along 
the line of sight. Cosmic defects, as well as open uni- 
verse or A dominated models provide explicit examples 
of theories where this actually happens. 

The smaller angle anisotropies are a more promising 
probe because they are due to local effects on the sur- 
face of last scattering, which are strongly constrained 
by causality (Figure 1). In particular, the pattern of 
Doppler peaks caused by phase coherent oscillations 0] 



FIG. I. The causality constraint on the microwave back- 
ground anisotropy. The picture is in comoving coordinates: 
the outer circle represents our causal horizon in the standard 
big bang. The inner circle is the surface of last scattering, 
on which photons are 'set free' from the hot plasma on their 
path to us. Circles show the domains of influence on these 
photons - the radius is the light travel distance since the hot 
big bang tls- This subtends an angle Ols ~ 1-1° in a flat 
universe with standard recombination. No causal physics op- 
erating within the standard big bang could have generated 
correlations between photons on the last scattering surface at 
points separated by more than 2€>ls on the sky. 



in the photon-baryon fluid provides a possible signature 
of 'super-horizon' curvature perturbations. Early papers 
exploring this include refs. j|, (§ and §. This Letter 
is a follow-up of U , where an inflation- mimicking model 
was proposed. A recent paper by Hu, Spergel and White 
criticising that model is now undergoing revision. 

In this Letter I ask whether a causal source acting 
purely via gravity within the standard hot big bang could 
generate CMB anisotropies similar to those in flat infla- 
tionary models. I further restrict consideration to 'scal- 
ing' sources, mainly for simplicity. This restricted set 
of models turns out to provide a surprisingly simple in- 
flationary mimic. I emphasise that the mimic is not a 
theory, but simply an 'ansatz' constructed by hand to 
provide a consistent solution to the Einstein equations. 
As such it provides a counterexample to some arguments 
in the literature (see e.g. Q). Furthermore it has a suf- 
ficiently simple form that it is not out of the question 
that it might actually be realised in a future theory of 
structure formation - in that sense the counterexample 
may turn out to be constructive. 
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I deal with the linearised Einstein equations in the 
'stiff' approximation, in which the induced perturbations 
are assumed to have negligible effect on the source . In 
this approximation, the source stress energy tensor Mt/ 
is covariantly conserved with respect to the background 
metric: 

e 00 + -(e 00 + e) = n ; ri + 2-n = c^e y (i) 

a a 

where IT = diQoi- Dots denote derivatives with respect 
to conformal time t, and a(r) is the scale factor. 

A formalism for dealing with such sources was pro- 
posed in Q. Here I shall consider only 'coherent' sources, 
which are representable in terms of a single set of 'mas- 
ter functions'. These are ordinary functions, which obey 
(|l|). The unequal time correlator < 6 Ml/ (r, t)9 p a(0, r') > 
equals the spatial convolution of the master function for 
Mt/ with that for <d p \. 

As argued in the causality constraint implies that 
the master functions Q^ v {r, r) are zero for all r > r. If we 
ignore vector and tensor perturbations, they are spheri- 
cally symmetric, and we can represent them as &oo(r, r), 
©oi = XiJ(r,r), Qij = ^0(r, r)<% + Q A (r,T)(x i x J - 
\r 2 5ij). The term ^0 is the pressure P, and <d A the 
anisotropic stress. In situations where matter is being 
actively moved, as it will be here, the anisotropic stresses 
are generally of the same order as the pressure. 

Some general properties can now be seen. The Fourier 
transforms (assumed to exist) are analytic about fc; = 0, 
and can be expanded in a Taylor series. The leading 
terms are then (by isotropy) 9oo ~ k , ©oi ~ h (from 
which n ~ k 2 follows) and 0^- ~ 8{j. In Fourier space 
we write 0ij(k) = |<%0 + (kjuj — ^Sij)Q s , where s = 
k- 1 dB A /dk - d 2 Q A /dk 2 . It follows that 5 - k 2 . 

I now specialise to 'scaling' sources, in which one as- 
sumes that the source a) involves a number with dimen- 
sions of the inverse of Newton's constant G and b) in- 
volves no other length scale apart from the horizon scale 
r. If these conditions are met, the source-perturbation 
equations are scale- invariant, apart from the violation of 
scaling caused by the radiation-matter transition. Scal- 
ing and dimensional analysis imply that (see e.g. j^]) 
0oo(fc,r) ~ r-i/i(fcT), 0(fc,r) ~ 7-i/ a (fcr), H(k,r) ~ 
fa(kr), s (fc,r) ~ T~2fa(kr), where the fa have 
power series expansions in k 2 obeying the restrictions 
noted above. These four fas are are related by the two 
energy momentum conservation equations (0). So for 
example, the k° term in the first equation relates the 
leading terms in fa and fa, and the k 2 term in the sec- 
ond equation relates the leading terms in fa and fa. But 
even after applying these equations, we still have essen- 
tially two free functions remaining. We also have some 
freedom in how to incorporate the matter-radiation tran- 
sition into the source. 

I assume the background spacetime is flat, and has 
metric ds 2 = a 2 (r)(— e?r 2 + (Sij + hij(x,T))dx i dx : >), with 



r conformal time and a(r) the scale factor. I work in ini- 
tially unperturbed synchronous gauge, which is especially 
suitable for causal theories since the Einstein equations 
are manifestly causal (in contrast, in Newtonian gauge 
the influence of anisotropic stresses propagates acausally 
||). This gauge is also straightforward to interpret - for 
example the radiation density contrast Sr is just that 
measured in the rest frame of freely falling particles like 
cold dark matter particles. 

We are interested in computing the CMB temperature 
distortion on the sky in a direction n: 

ST 1 1 f f 

— (n) = -5 R {i) - n • v R (») - - J drh^nV (2) 

where Sr is the density contrast, and Vr the velocity, of 
the photon fluid on the surface of last scattering. The 
last term is the Sachs- Wolfe integral, representing the 
change in the proper path length field along the line of 
sight. This expression holds in the 'instantaneous recom- 
bination' approximation, which I shall use throughout. 

The first two terms are local effects on the surface of 
last scattering. They are determined from: 

S'c + -5c = 4^G(V(1 + 3c 2 N ) PN S N + 0oo + ©), (3) 

N 

Sr = ^S c - • v R ;v R = -(1 - 3c|)^v R - ^c 2 s X7S R 

(4) 

where cs is the speed of sound in the photon-baryon fluid. 
The only component of the stress energy tensor that en- 
ters is ©oo + ©, i-e. p + 3P. Thus prior to last scattering, 
only one of the two free functions in 0^^ contributes - 
the other is literally 'invisible' in the CMB anisotropy. 

In the simplest inflationary theory, the surface term 
jSr captures most of the relevant physics determining 
the location of the Doppler peaks. In ref. Q, I focussed 
only on this term, and found a causal source 0oo + © 
which mimicked the same term in the inflationary theory. 
The idea was simply to choose 

eoo + e<x/i(r)+/ a (r)«*(r-i4T) 0<A<1 (5) 

representing a spherical shell expanding at some frac- 
tion of the speed of light. I found that for A close to 
unity, the Sr surface term closely matched the inflation- 
ary one. Such a shell of matter is similar in form to a 
supernova explosion - for a spherical shell of neutrinos, 
one has ©^ ~ Y,-p % p> oc x l x^ jr 2 . In that case, s and 
are comparable in magnitude, and the same is true here. 

Here I extend the computation to the entire expression 
(^). The Sachs Wolfe integral introduces some depen- 
dance on the anisotropic part of the metric perturbation. 
To compute this it is necessary to further specify 0/^, 
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by for example sepcifying another of the fa functions. 
The simplest choice leaving /i + / 2 fixed is to specify fz- 
Then equations (|]) are used as follows: the energy equa- 
tion is integrated to determine 9oo, and the momentum 
equation is differentiated to determine 9 5 . Of course 
this must be done consistently with the matching of the 
leading terms as discussed above. 

In Fourier space, the choices I make for the source are: 
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oo 



a sim4fcr 
a Akri 



(0) 



as in (|^) , with the prefactor incorporating the radiation- 
matter transition in a simple way. For II, we must sat- 
isfy II(fc) ~ k 2 at small k. Equivalently, the integral 
J r 2 drH(r,T) = 0. This is most easily satisfied by tak- 
ing II(r, t) to be the sum of two delta functions, of equal 
weight but opposite sign. Their Fourier transform pro- 
duces: 



n 



E(t) 6 ,sinBkT smCkr^ 



kri B 2 - C 2 



D 



C 



(7) 



where E(t) is a messy function obtained by analytically 
solving for the coefficient of fc 2 term in the momentum 
equation (Jl]). It equals j§ in the radiation era and jg 
in the matter era. A set of values which I find leaves 
the Sachs Wolfe integral sub-dominant is B = 1.0 and 
C = 0.5. 

The anisotropic metric perturbation h s (k) (defined 
analogously to 9 s ) is now given by h s — h = — 247rG(n + 
^n(Pn + /5Ar)a 2 ik • vj\r)/fc 2 . Finally, I model the free 
streaming of photons and neutrinos after last scattering 
following ref. ||. The initial conditions for the 9oo and 
the perturbations are set up deep in the radiation era 
t << teq, and well outside the horizon fcr << 1: they 
are read off from the k° terms in the energy conservation, 
and perturbation equations: 



9c 



2r" 



,8c 



v fi = (8) 



with the constant D determined by setting the to- 
tal pseudoenergy t o = k 2 {h — /i s )/(247rG) = 9oo + 
^2 N pno 2 8m+{o,/ a)5c I (4:irG) to zero. With these choices 
there are no superhorizon perturbations in the photon- 
to-CDM, photon-to-baryon or photon-to-neutrino ratios. 
The pseudoenergy is just the Ricci scalar of the spatial 
slices, so setting it zero means there are no curvature per- 
turbations either. These initial conditions are thus both 
'adiabatic' and 'isocurvature' - there are simply no per- 
turbations in the universe on superhorizon scales. The 
complete Ci spectrum of the causal model defined in 
equations (6-8) is shown in Figure 2. 

What about the sub-horizon behaviour of the source? 
In the construction above, where I integrate the energy 
equation to determine 9oo, there is no reason it should 
tend to zero well inside the horizon. However, because I 
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10 100 1 1000 

FIG. 2. Comparison of the simplest inflationary theory 
with its 'mimic' causal source model discussed here. The 
vertical axis is 1(1 + l)Ci, with Ci the angular power spec- 
trum, and I the Legendre index. Both curves were calcu- 
lated in the same 'instantaneous recombination' approxima- 
tion, in a flat universe with 'canonical' parameters SIb = 0.05, 
FIcdm = 0.95. The vertical scale is arbitrary. 



explicitly turn off 9oo + 9 inside the horizon, the source 
ceases to have any effect on the fluid perturbations and 
the trace part of the metric h = —25c- The effect on 
the anisotropic part h s is similarly turned off because II 
goes to zero. In effect, I have turned off all the 'grav- 
itationally active' components of the source, but there 
is no reason for the energy 9oo, the pressure 9 or the 
anisotropic stresses 9 s to separately vanish - they only 
need to satisfy the relations 9oo+9 = and 9+29 s = 0. 
This is reminiscent of the behaviour of a straight cosmic 
string - it carries energy but generates no gravitational 
field. In any case, the sub-horizon source is removed if 
one adds a term c\k 2 T(a/a)Q^ to 9oo + 9, or a term 
— C2& 2 t9oo to II. Either of these makes 9oo, 9, and 9 s 
go to zero inside the horizon as exp(— const. k 2 T 2 ). Fig- 
ure 3 shows the evolution of the components of with 
and without these modifications, and Figure 4 shows the 
corresponding Ci spectra. The moral is that there is a 
lot of freedom inside the horizon to make very large al- 
terations in the source without significantly affecting the 
CMB anisotropies. 

How sensitive to the particular choice of Ansatz is the 
result? Figure 4 compares the ^4 = 1 model with the 
cases A — 0.7, and A — 0.1. The leftward shift in the 
first peak was explained in ref. |@[. As shown there, it 
is also easy to arrange for a shift to the right. So while 
the A = 1 model is on the causality limit, there is a sub- 
stantial region of parameter space around it with similar 
Ci spectra. It would be very difficult to observationally 
separate these models from inflation, especially in the 
realistic case where we are unlikely to know all the rele- 
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FIG. 3. Behaviour of the stress tensor inside the horizon. 
The thin lines show the original source denned in equations 
(6) and (7), as a function of kr in the radiation era (in order 
as the curves intersect the y axis moving upwards) r 5 0s, 
7-i(0 + 0oo), 7-i0oo and 75k~ 2 r~ill. The bold lines show 
the same curves for the model specified by ci = .001 and 
C2 = .003 ( see text), where all components are 'turned off' 
inside the horizon. The Ci spectrum for the latter model is 
shown in Figure 4. 
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FIG. 4. The Ci spectrum of the A = 1 model (bold line) 
is compared with that for A = 0.7 (dashed line) and A = 0.1 
(dotted line). The model illustrated in Figure 3, for which 
the stress tensor vanishes inside the horizon, is also shown 
as the thin line (altering C2 to 0.001 makes the Ci spectrum 
virtually indistinguishable from the original model, except for 
the fourth peak, which is still a little high)). 



vant cosmological parameters (fi, fi^, h, fl\ and so on) 
in advance. I have so far explored only a very tiny part 
of model space, and preliminary investigations of more 
general Ansatzes indicate that a very wide range of Ci 
spectra are possible, especially when the Sachs Wolfe in- 
tegral becomes a dominant effect. 

I have given myself a great deal of freedom in con- 
structing this source, and it is far from clear that realistic 
physics could produce it. Nevertheless, the reader may 
tolerate some speculations. As mentioned, the form of 
p + 3P required is similar to that resulting from a super- 
nova explosion, but whereas in the latter case the energy 
redshifts away as a(t) -1 , here, scaling evolution requires 
that the energy in the shell increases - by dimensions, 
E ~ G~ 1 t. This requires some positive feedback mecha- 
nism, which one could conceivably arrange with unstable 
dark matter, decaying via stimulated emission of Gold- 
stone bosons, or even gravity waves. Another issue is 
the Gaussianity or otherwise of the perturbations. If the 
source is made up of a very large number of 'explosions' 
which arc allowed to superpose, then it can be made ar- 
bitrarily Gaussian. But if the exploding shells interact, 
there would be a limit to their number density, and the 
perturbations would then be nonGaussian. 

The conclusion of this paper is that causality alone is 
insufficient to distinguish the inflationary Ci predictions 
from those of non-inflationary models. Of course the ob- 
servational confirmation of one of these spectra would be 
a tremendous success for inflation, but the door would 
nevertheless still be left open to other possible explana- 



tions of cosmic structure formation. 
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